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In a simultaneous paired approach to closely-spaced parallel runways, a pair of 
aircraft flies in close proximity on parallel approach paths. The longitudinal 
separation between the aircraft must be maintained within a range that avoids wake 
encounters and, if one of the aircraft blunders, avoids collision. To increase 
operational availability, the approach procedure must accommodate a mixture of 
aircraft sizes and, consequently, approach speeds. In these procedures, the slower 
aircraft is placed in the lead position. The faster aircraft maintains separation from 
the slow aircraft in a dependent operation until final approach and flies 
independently afterward. Due to the higher approach speed of the fast aircraft, 
longitudinal separation will decrease during final approach. Therefore, the fast 
aircraft must position itself before the final approach so that it will remain within 
the safe range of separation as separation decreases. Given the approach geometry 
and speed schedule for each aircraft, one can use kinematics to estimate the 
separation loss between a pair of aircraft. A kinematic model can complement fast- 
time Monte-Carlo simulations of the approach by enabling a tailored reduction in 
the variation of starting position for the fast aircraft. One could also implement the 
kinematic model in ground-based or on-board decision support tools to compute the 
optimal initial separation for a given pair of aircraft. To better match the auto- 
coupled flight of real aircraft, the paper derives a kinematic model where the speed 
schedule is flown using equivalent airspeed. The predicted time of flight using the 
equivalent airspeed kinematic model compares well against a high-fidelity aircraft 
simulation performing the same approach. This model also demonstrates a modest 
increase in the predicted loss of separation when contrasted against a kinematic 
model that assumes the scheduled speed is true airspeed. 
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Nomenclature 

glideslope angle (radians) 

atmospheric density at altitude (slug per cubic foot) 
atmospheric density at sea-level (slug per cubic foot) 

Automatic Dependent Surveillance-Broadcast 

constant acceleration of equivalent airspeed, negative on approach (knot/s) 
above ground level 

acceleration of the slow aircraft, negative on approach (knot/s) 

Closely Spaced Parallel Operations 

compression of longitudinal separation as the slow aircraft travels from FAF to PTS 
collision safe distance (feet) 
equivalent airspeed distance (feet) 

true or D EAS distance traveled by the aircraft during deceleration (feet) 

Deas of the fast aircraft from PTF when the slow aircraft is at the FAF and deceleration is 
independent (feet) 

Deas of the fast aircraft from PTF when the slow aircraft is at the FAF and deceleration is 
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dependent (feet) 

offset distance between runways (feet) 
wake-safe distance (feet) 
equivalent airspeed (knots) 

conversion factor for EAS to TAS as a function of geometric altitude 
final approach fix 
feet per second 

standard acceleration due to free fall (ft/s 2 ) 

geometric altitude, assumed same as altitude above mean sea level (feet) 
equivalent airspeed altitude (feet) 
geopotential height (feet) 

geometric altitude of the final approach fix (feet) 

geometric altitude of the fast aircraft when the slow aircraft is at the FAF (feet) 

geometric altitude of the procedure termination point for the fast aircraft (feet) 

geometric altitude of the procedure termination point for the slow aircraft (feet) 

geometric altitude of the stabilized approach point (feet) 

geometric altitude of the threshold crossing height (feet) 

initial approach fix 

knots 

temperature lapse rate (Rankine per foot ) 8 

atmospheric pressure at altitude (pound force per square foot) 

atmospheric pressure at sea level (pound force per square foot) 

pound force per square foot 

procedure termination point of the fast aircraft 

procedure termination point of the slow aircraft 

gas constant for the atmosphere 

equatorial radius of the Earth (feet) 

stabilized approach point 

required longitudinal separation at the FAF (ft) 

atmospheric temperature at altitude (Rankine) 

atmospheric temperature at sea-level (Rankine) 

true airspeed 

threshold crossing height (feet) 

time for an aircraft to decelerate from V c to V f (seconds) 
time delay in the dependent operation response of the fast aircraft (seconds) 
flight time from the FAF to the SAP (seconds) 
flight time of an aircraft in the final approach segment (seconds) 
flight time of the fast aircraft from FAF to PTF under independent operation (seconds) 
flight time from the SAP to the PTF (seconds) 
flight time from the SAP to the PTS (seconds) 
flight time from the SAP to the TCH (seconds) 
flight time from the FAF to the PTS for the slow aircraft (seconds) 
equivalent airspeed at start of deceleration (knots) 
constant segment speed (knots) 
equivalent airspeed (knots) 
final approach speed (knots) 
true airspeed (knots) 

wind speed as a function of altitude (knots) 
the x-axis runway location of the fast aircraft when the 
fast aircraft deceleration is dependent (feet) 
the x-axis runway coordinate of the FAF (feet) 
the x-axis runway location of the fast aircraft when the 
fast aircraft deceleration is independent (feet) 


slow aircraft is at the FAF and the 


slow aircraft is at the FAF and the 
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I. Introduction 

T HE Closely Spaced Parallel Operations (CSPO) Working Group of the Federal Aviation 
Administration is investigating operational concepts to allow simultaneous paired approaches under 
instrument conditions at runways separated by as little as 700 ft. 1,2 In a simultaneous paired approach, a 
pair of aircraft flies in close proximity to each other as they approach a pair of parallel runways. The 
longitudinal separation between the aircraft must be maintained within a range to avoid wake encounters 
and, if one of the aircraft blunders^, to avoid collision. To increase operational availability, the approach 
procedure must accommodate a mixture of aircraft sizes and, consequently, approach speeds. The aircraft 
with the slower approach speed is placed in the lead position and the faster aircraft in the trail position. The 
fast aircraft first establishes dependent operation to the slow aircraft. The fast aircraft maintains speed and 
position against the slow aircraft until the fast aircraft reaches its final approach speed, which occurs after 
passing the final approach fix (FAF). For the remainder of the approach, the faster aircraft flies 
independently. It will close in on and, if the procedure allows, possibly pass the slower aircraft before 
landing. To appropriately position the fast aircraft before the FAF, it is necessary to determine the amount 
of separation that will be lost during the final approach. 

The approach speeds of each aircraft and the approach geometry are the primary determinants of the 
separation loss during the approach. Given the approach geometry and speed schedule for each aircraft, 
one can use kinematics to compute the separation loss between a pair of aircraft. A kinematic model can 
complement Monte-Carlo simulations of the approach by enabling a tailored reduction in the variation of 
starting position for a given aircraft pair through identification of the ideal initial separation. The kinematic 
model could also be inserted into ground-based or on-board decision support tools to compute the optimal 
initial separation for a given pair of aircraft. 

The precision of the kinematic model is limited by its assumptions. One challenge to improving the 
precision of a kinematic model is that aircraft do not fly the speed schedule using true airspeed. Aircraft 
use indicated or calibrated airspeed, which are functions of atmospheric properties at altitude. Therefore, 
even in segments of the descent where the scheduled speed is constant, the aircraft executes a slow 
deceleration during descent. Because the initial longitudinal separation of aircraft in a paired approach is 
less than a nautical mile, the impact of the slow deceleration on loss of separation, though small, is 
significant. Therefore, the kinematic model should capture this effect. Using the 1976 U.S. Atmosphere 
Model, the paper formulates a kinematic model using equivalent airspeed (EAS) to better approximate the 
indicated or calibrate airspeed flown by the aircraft. That EAS model mirrors a kinematic model that uses 
true airspeed by introducing an equivalent airspeed distance (D EAS ). The predicted loss of separation using 
the EAS model is compared against a high-fidelity simulation performing the same approach. 

II. Approach Geometry 

Candidate geometries for closely spaced parallel approaches (CSPO) are detailed in references 1 
through 6. To summarize, air traffic control vectors a pair of aircraft onto parallel approach paths at a point 
between the initial approach fix (IAF) and the final approach fix (FAF). The aircraft with the faster final 
approach speed is placed behind the aircraft with the slower final approach speed. The aircraft also begin 
with a vertical separation of 1000 ft to allow each aircraft to safely establish their approach path and to 
allow the fast aircraft to establish longitudinal separation from the slow aircraft. The faster aircraft can 
start at either the higher altitude or the lower altitude. Each aircraft flies at a constant altitude until 
encountering the glidepath and then descends along the glidepath. From the start of the procedure to the 
FAF, the slow aircraft flies an assigned constant speed (nominally 180 KT) and the fast aircraft maintains 
speed and position against the slow aircraft. The fast aircraft executes this dependent operation until it 
reaches its final approach speed (matching speed of the slow aircraft as the slow aircraft decelerates upon 
reaching the FAF) or, if the fast aircraft does not detect the acceleration of the slow aircraft prior to the 
reaching the FAF, the fast aircraft switches to independent operation at the FAF. 


^ An aircraft is said to blunder when it departs from its approach path toward the path of the other aircraft. 
A blunder may occur due to pilot error or aircraft malfunction. 
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Figure 1: Paired Approach Separation Distances 

Figure 1 depicts the horizontal geometry of the paired-approach procedure. The figure shows two 
options for the separation window (i.e. the green box), the full length of the box applies to a procedure that 
allows the fast aircraft to pass the slow aircraft 1 and the hashed boundary applies to a procedure that does 
not allow the fast aircraft to pass 6 . The former requires the trail aircraft to be equipped with a collision 
alerting and avoidance system that is designed for the tight vertical and longitudinal separations in the 
paired approach. The latter is designed to be passively safe against collisions and, therefore, requires no 
special equipage. In each option, the separation window is defined to avoid two hazards. In the case where 
passing is not allowed, those hazards are collision and wake. To avoid collision, the fast aircraft must 
maintain a sufficient distance behind the slow aircraft such that the slow aircraft will not collide with the 
fast aircraft if the slow aircraft blunders across the fast aircraft's path. To avoid wake, the fast aircraft must 
be close enough to the slow aircraft that the wake of the slow aircraft passes behind the fast aircraft in the 
presence of adverse crosswind. In the case were passing is allowed, the two hazards are both wake 
encounters. The fast aircraft must remain forward of a rear wake-safe boundary to avoid encountering the 
wake of the slower aircraft. The fast aircraft must also remain behind a forward wake -safe boundary so 
that the slow aircraft does not encounter the wake of the fast aircraft. The paired approach with passing 
relies on collision alerting rather than separation to prevent collisions when one or both aircraft blender. 
With or without passing, the fast aircraft must maintain its separation between the forward and rear 
boundaries from the loss of altitude separation until one of the aircraft crosses the threshold. ^ 

The kinematic models derived in this paper assume that the runway thresholds are not offset but they 
can have different elevations. Section III.E.3 discusses adaptation of the models for offset runways. If the 
runway thresholds are offset, then the slow aircraft is assigned to the closer runway. From inspection, the 
kinematic models should be applicable to runways with different glideslope angles, but the models have not 
been verified and validated under this scenario. If a runway pair does have differential glideslopes, then the 
slow aircraft flies the shallower glideslope since wake vortices tend to descend. 


* Both collision and wake risks are assumed to be negligible after the lead aircraft crosses the threshold. 
The assumption on collision risk is reasonable given that ground hazards restrict banking by the lead 
aircraft and the go-around procedure can be designed to increase lateral separation. However, further wake 
studies on landing and go around scenarios are needed to validate the assumption of negligible wake risk. 
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III. Modeling Travel Time and Separation Loss 

Prior to the FAF, the fast aircraft uses speed management to actively control its longitudinal separation 
with the slow aircraft. The fast aircraft stops active management of longitudinal separation when it reaches 
its final approach speed (by decelerating in concert with the slow aircraft after the slow aircraft reaches the 
FAF) or when it reaches the FAF (if it has not detected deceleration by the slow aircraft beforehand). For 
the remainder of the procedure, the longitudinal separation will decrease due to the difference in approach 
speeds. However, the fast aircraft continues to monitor the longitudinal separation and will execute a go- 
around if it falls outside the separation window. Therefore, to limit nuisance go-arounds, the fast aircraft 
must position itself to accommodate the longitudinal compression before the slow aircraft reaches the FAF. 

If the aircraft accurately fly their planned approach speed to the threshold, then the longitudinal 
compression is characterized by the approach speeds of the two aircraft, the rate of deceleration, the 
minimum separation that must be maintained, the height to which minimum separation must be maintained, 
winds, and atmospheric density. Winds and atmospheric density come into play because the approach 
speed is an indicated (or calibrated) airspeed, not a true ground speed. For simplicity, the model uses 
equivalent airspeed to approximate indicated airspeed. Also, winds will be ignored for the moment and 
addressed later in the paper. 

Accounting for atmospheric density with altitude is desirable. Assuming a 1976 standard atmosphere, 
the equivalent airspeed at an altitude of 1800 ft is about 3% lower than the true ground speed. From a FAF 
at 1800 ft altitude to a touchdown crossing height (TCH) at 50 ft atltitude, the elapsed time to threshold is 
about 2 seconds shorter for the scheduled speed as EAS than if one assumes the scheduled speed is TAS. 
The challenge, however, in computing compression distance using EAS is that the conversion of EAS to 
TAS (necessary to apply kinematic principles) changes with altitude as density changes with altitude. 

A. Model with TAS = EAS Assumption (TAS Model) 

Before incorporating the EAS to TAS conversion, the kinematic equations for the approach procedures 
are first formulated assuming the speed schedule is TAS (TAS model). Later, this paper will show how to 
formulate the EAS to TAS conversion in a way that resembles the TAS model. The speed schedule for the 
paired approach is depicted in Figure 1. The slow aircraft's speed schedule can be broken into three 
segments: a constant speed segment, a deceleration segment, and a final approach segment. The constant 
speed segment begins when the slow aircraft captures the glideslope and ends at the final approach fix. In 
this region, the aircraft flies a constant airspeed assigned by air traffic control. When the slow aircraft 
reaches the FAF, it begins the deceleration segment where it decelerates to its final approach speed. The 
model assumes a constant deceleration that achieves the final approach speed at the stabilized approach 
point (SAP), i.e. 1000 ft AGL. The final approach segment continues after the SAP; the slow aircraft flies 
its final approach speed to the runway threshold. Separation monitoring ends when one of the two aircraft 
reaches the threshold so there is no need to model the flare maneuver after crossing the threshold. 

To compute the compression distance 
at the FAF, one must first calculate the 
time it takes the slow aircraft to travel from 
the FAF to the end of the paired approach. 

If the fast aircraft is not allowed to pass, 
then the paired approach ends for the slow 
aircraft at the threshold. If the fast aircraft 
is allowed to pass, then the paired approach 
ends when the fast aircraft crosses the 
threshold. To avoid the wake of the fast 
aircraft, the slow aircraft must be behind by 
no more than the wake-safe distance (D ws ). 

The location of the slow aircraft at the end 
of the paired approach will be called the 
procedure termination location for the slow 
aircraft or PTS. 

The flight time from FAF to PTS is 
divided into two computations: the time 
from FAF to SAP and the time from SAP 


Distance From Threshold 

► 

FAF 



Figure 1. Speed Profiles of Slow and Fast Aircraft 
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to PTS. The time from FAF to SAP assumes a constant deceleration from the constant segment speed (V c ) 
to the final approach speed (V f ). Thus, the average speed from FAF to SAP is V 2 [V c + V f ]. Given that the 
FAF and SAP are defined by height above the runway, the distance between these points is a function of 
their altitudes and the glideslope angle: 


= DQi,y) = (hfAF - hsAp) 2 KT = 18113 
Vavg siny (V c + V f )fps (V c + V f ) 


( 1 ) 


where, 


D(h,y) is the distance traveled as a function of altitude and glideslope angle 

V avg is the average true speed of the aircraft 

h FAF is the geometric altitude of the FAF in feet, normally 1800 ft above runway elevation 

h S AP is the geometric altitude of the SAP in feet, normally 1000 ft above runway elevation 

y is the glideslope angle, normally 3° 

V c is the constant segment speed in knots, nominally 180 KT 

V f is the final approach speed in knots 

KT/fps is the conversion from knots to feet per second, 0.5924838 
t F AF->SAP is the flight time from the FAF to the SAP in seconds 


Similarly, the time from SAP to PTS is given by: 


_ DQi,y) _ ( h S AP ~ h PTS ) KT 
Vf Vf sin y fps 


( 2 ) 


where, h PTS is the altitude of the procedure termination location for the slow aircraft: 

hpTs = hpcH + A vs t an Y (3) 


InEq. (3), 

D ws is the wake-safe distance. If passing is not allowed, this variable is zero. 

h T cH is the altitude of threshold crossing height, nominally 50 ft above runway elevation 


In the example of a paired approach without passing, D ws is zero and h PTS equals h TCH . Given a generic 
runway at sea-level elevation with a TCH of 50 ft AGL, Eq. (2) reduces to: 


= ( hsAP ~ h TCH ) KT = 10755 
Vf sin y fps Vf 


(4) 


Then, the total flight time from FAF to TCH for the slow aircraft landing at the generic runway is: 


t 


slow 


t FAF -> SAP + tsAP^PTS 


18113 ( 10755 

( Vc + Vf) + Vf 


(5) 


For example, if V c is 180 KT and V f is 120 KT, then the t s i ow is 150.0 seconds. 

Next, one calculates where the fast aircraft would be t s i ow seconds in the past. This calculation depends 
on the separation that the fast aircraft must maintain relative to the slow aircraft at the runway threshold. In 
the paired procedure without passing, the faster aircraft must maintain a collision free distance (D C s) from 
slower aircraft upon procedure termination. For the paired approach with passing, the procedure terminates 
for the fast aircraft at the threshold. As with the slow aircraft, this required separation is converted into a 
procedure termination altitude for the faster aircraft, h PTF : 

hpTp = h TCH + D cs tan y (6) 


where, 

D cs is the collision-safe distance. If passing is allowed, this variable is zero. 

h T cH is the altitude of the threshold crossing height, nominally 50 ft above runway elevation 

y is the glideslope angle, nominally 3° 
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Note that the h TCH in Eq. (6) can differ from h TCH in Eq. (3) to accommodate runways with different 
threshold crossing heights; moreover, because the geometric altitude and not the runway relative altitude is 
used, the runways can also have different elevations. Similarly, for Eq. (7) to Eq. (15) that follow, h S AP and 
h FAF can differ for the fast aircraft to accommodate different runway elevations or approach geometries. To 
handle differential glideslopes, the glideslope angle in Eq. (7) to Eq. (15) can also differ from that in Eq. (3) 
except where noted. 

One aspect of the paired procedure that complicates the computation is the condition on the time and 
deceleration for the fast aircraft. The fast aircraft decelerates at the same rate as the slow aircraft when and 
if it detects the deceleration of the slow aircraft prior to the fast aircraft reaching the FAF; otherwise, the 
fast aircraft decelerates independently. The parameter governing this condition is t de i ay , the time delay in 
the response of the fast aircraft to changes in the sensed state of the slow aircraft. This exercise will 
assume that the speed of the slow aircraft is monitored using ADS-B and exhibits a t de i ay of 5.0 seconds. 
Table 1 provides a breakdown of the estimated t de i ay . 


Table 1. Estimate of Response Delay in Dependent Operations 


ADS-B OUT Fatency (from 
navigation measurement) 7 

2.0 s 

ADS-B IN Fatency 

0.5 s 

Speed Guidance Algorithm 

0.5 s 

Pilot Reaction Time 

2.0 s 

Total 

5.0 s 


Let q.fast be the independent flight time of the fast aircraft from the FAF to PTF. The fast aircraft will 
decelerate independently if (t sk)W - h.^) < t de i ay . The equations for t T . fast are identical to those for t s i ow but 
h PTF is substituted for h PTS . For this example using the generic runway, the collision-safe distance, D C s, is 
set to 750 ft. 

h PTF = h TCH + D cs (tany ) = 50 ft + 750 ft (tan 3°) = 89.3 ft (7) 


_ D(h,y ) _ (h SAP - h PTF ) KT _ 10310 
V VfSiny fps V f 


( 8 ) 


18113 ( 10310 

tl—fast ~ ^FAF^SAP + ^SAP^PTF ~ r, T " \ 3 7 } (9) 

(Vc + Vf) Vf 


If the value of ti_ fSLSt reveals that the fast aircraft will perform an independent deceleration, then the position 
of the x-axis runway coordinate of the fast aircraft when the slow is at the FAF is: 


X i- fast 


= Xfaf - cosy 




( 10 ) 


where X FAF is the x-axis runway coordinate of the FAF, nominally -33,392 ft or -5.5 NM. If X FAF differs 
between the two runways, then use X FAF for the runway of the fast aircraft in Eq. (10). 

If the fast aircraft will perform a dependent deceleration with the slow aircraft, the fast aircraft's total 
flight time, t D . fast , can be broken into three components: 1) t de i ay - the delay in dependent operation defined 
in the previous paragraph, 2) t dece i - the time to decelerate to the final approach speed, and 3) t fina i - the time 
flying at the final approach speed. t dece i is computed using the deceleration of the slow aircraft (due to the 
dependent operation) and the final approach speed of the fast aircraft: 


_ / \ - iu v ■) + V A sin 7 ill ~ z V H 

slow \t p af^ S ap si 0 w f c 2 (h FAF — h SAP ) KT 18113 


( 11 ) 
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fast 


( 12 ) 


slow 


Note that the equation for a s i ow produces deceleration in units of kt/s so that the equation for t dece i produces a 
time in units of seconds. Also, note that in the case of differential glideslopes, Eq. (11) uses the glideslope 
angle for the slow aircraft. Now, t fina i becomes the remaining time: 


t final t 


slow t delay t, 


decel 


Next, compute the runway x-axis location of the fast aircraft performing a dependent deceleration: 
X D-fast ~ ~ ^CS + COS y {Vf tfinai ^/) tdecei Vc tdelay ) ) 

Now, compute the required longitudinal separation at the FAF for the aircraft pair: 

i % I-fast if ( tslow ~ ll-fast) ^ ^delay 


$FAF — XpAF 


X D-fast if ft s i ow li-fast) — ^delay 


(13) 


(14) 


(15) 


If X FAF differs between runways, then use X FAF for the runway of the slow aircraft in Eq. (15). Equation 
(15) assumes that the runway thresholds are not offset in the x-axis direction. 

The compression distance is the difference between longitudinal separation at the FAF and longitudinal 
separation at the threshold: 


* compress 


— $faf — 


cs 


+ A 


ws 


(16) 


(Recall that D cs is zero for the paired approach with passing and that D ws is zero for the paired approach 
without passing. So, at procedure termination, the fast aircraft is either a distance of D cs behind the slow 
aircraft or a distance of D ws ahead of the slow aircraft.) 

Continuing the example, if one computes the compression distance for a paired approach without 
passing where the slow aircraft final approach speed is 120 KT, the fast aircraft final approach speed is 130 
KT, the constant segment speed is 180 KT, and the collision safe distance D C s is 750 ft, then the results are: 
t s i ow = 150.0 seconds 

h-fast =137.7 seconds, (t sk)w - ti_ fast ) =12.3 seconds > t dday so compute X D _ fast 

t d ecei = 50.3 seconds 

X D _ fast = -36159 ft 

S faf = 2767 ft 

^compress — 2017 ft 

B. Including the EAS to TAS Conversion (EAS Model) 

Next, one incorporates the EAS to TAS conversion to fly the speed schedule as EAS: 



where p, P, and T are the atmospheric density, pressure, and temperature respectively and the subscript zero 
indicates the sea-level value. The second formulation allows the insertion of functions for temperature and 
pressure from the US Standard 1976 atmosphere model for altitudes below 36,152 ft: 8 


H 


hR e 
h + R e 


(18) 
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(19) 


T = 


T 0 (l 


H L\ 

+ T> 


P 


/ HL\ 

p °( 1+ id 


-g/(L R ) 


( 20 ) 


where, 

g = the standard acceleration due to free fall, 32.17404 ft/s 2 

h = the geometric altitude, feet. Equal to mean sea level altitude. 

L = the lapse rate, -3.56616xl0‘ 3 °R/ft (converted from SI units in Ref. 8) 

P 0 = sea-level pressure, 21 16.22 psf 
R = the gas constant for air, 1716.55915670803 ft 2 /(s 2 °R) 

R e = the radius of the Earth, 20855531.5 ft 
T 0 = sea-level temperature, 518.67 °R 
H = the geopotential altitude, feet 


With these equations, the conversion from EAS to TAS becomes: 


VtAS ~ ^EAsfW 



hLR e \ 1+fl/(ifi) 
(h + R e )Tj 


= V, 


EAS 



143.394/1 > 

ft + 20855531.5/ 


- 4.25589 


( 21 ) 


Equation (21) computes TAS as a function of geometric altitude, h. To determine the distance traveled 
over a given time, we recognize that TAS is dr/dt where r is the distance along the glideslope. Since the 
glideslope angle is constant, we can convert dr/dt to dh/dt using the relation, h = r sin(y). 


-(■ 


/ KT \ dr / KT\dh 

‘ , ™ = fej* = lcscv — 


cscr TFs)^ = Vl ‘ sfm 


r /esc y KT W 1 

J dt 1 \Veas fps) J f(h) 


dh 


( 22 ) 

(23) 


The integral for dh does not have a closed form solution. However, one can construct a closed form 
approximation by converting the integrand to a Taylor series and integrating the first two terms: 


t = 



- 1.46309 x 10" 5 h + 5.74347 x 10" 11 ft 2 ) dh 


(24) 


cscy (ft - 7.31543 x 10" 6 ft 2 + 1.91449 x 10" 11 ft 3 ) KT D EAS (h) KT 


v EAS 


fps 


V, 


EAS 


fps 


As indicated by its final formulation, Eq. (24) defines an equivalent airspeed distance (D EAS ) that is a 
function of geometric altitude. Note that removing the second and third order terms reduces Eq. (24) to the 
assumption that TAS = EAS. Thus, the second and third order terms represent the effect of the EAS— >TAS 
conversion on the time to reach sea level. For example, when starting at a 1000 ft altitude, the higher order 
terms indicate that the EAS— >TAS conversion is equivalent to reducing the distance traveled by 139 ft. 
The third order term produces a distance of less than 1 foot below h = 3738 ft. It can be ignored for 
runways at elevations below 1500 ft. 
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We can substitute D EAS for the h/sin(y) terms in Eq. (1) through Eq. (12). For example, the time to 
travel from SAP to PTS, Eq. (1), becomes: 


l sap^pts 


DQhY) = (D EAS ( h S AP ) Peas (/i P75 )) KT 

Vf Vf 1 Vs 


(25) 


This substitution also works for the deceleration segment when deceleration is a constant deceleration 
of the equivalent airspeed because a constant deceleration term can be evaluated in the time integral: 

KT\dh 


Y 


TAS 


n ( KT\ Ctrl 

= (Yeas + “eas 0 f ( h ) = ( CSC ^^J ^ 


(26) 


f A fPS f 1 

J (Yeas + cLeas 0 ~£f dt = CSC T J J^ dh 


(27) 


(y EAst + 2 a EAS f 2 ) r^T — PeAsQ 0 


(28) 


Equation (28) is identical to the true airspeed kinematic equation for constant acceleration, with D EAS 
replacing the distance term. Finally, Eq. (1) through Eq. (12) are modified to use D EAS . Using example 
values for the generic runway, the final equations for a paired approach without passing are: 

. _ '(D eas (1000) - D eas (50)) 2(D eas (1800) - D EAS (1000))\ KT 

^slnw I , ^ ' - - - - I ‘ 


Vf 


f 


Vf + Vc 


J I fps 


10672 17744 

+ ■ 


V f 


v f + v c 


(29) 


f I- fast 


(D eas ( 1000) - D eas (50 + D cs tan 3°)) KT 17744 10228 17744 


fps + V f + V c 


■ + 


Vf + V c 


C^slow 


= {r^—) =( v f~ v c) 

^FAF^SAP' siow 


Vf + Y 

17744 


(Yf-Vc 2 ) 

17744 


(30) 

(31) 


AV fast 

taecel = ~r — = 17744 
Q-slow 


(Yf-Vc) 


fast 


(V/-Vc 2 ) 


slow 


(32) 


To compute the runway x-axis coordinate of the fast aircraft when the slow aircraft is at the FAF, one first 
computes the position as a D EAS distance: 


Pi -fast — PeAs(^faf) + Y(fslow 1 1 - fast ) 


fps 

KT 


(33) 


Pd -fast — DeaS QIptf) + tyf t final + ^(X + Y) t decel + Y t delay ) 

= 1705 + ( Vf tf ina i + ^(v c + Vf) t decei + V c t de i ay ) 


(34) 


D 


fast 


Pi-fast O' slow ll-fast) — ^delay 
Pp—fast if Oslow 1 1 -fast) ^ Idelay 


( 35 ) 


Then, solve the D EAS formula to determine the starting geometric altitude: 
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(36) 


esc/ (h - 7.31543 x 1(T 6 h 2 + 1.91449 x 1CT 11 h 3 ) = D fast 

1.91449 x KT 11 h 3 - 7.31543 x 1(T 6 h 2 + h- D fast sin / = 0 (37) 

The last term, D fast siny, represents an equivalent airspeed altitude, h'. Thus, the resulting solution will be 
a conversion from the equivalent airspeed altitude (h') to the geometric altitude (h). The real cube root of 
Eq. (37) is a complicated formula. To simplify it, the cube root is expanded into a Taylor series and the 
first four terms of the series are used: 

hfast = h' + 7.31543 x 1(T 6 h' 2 + 8.78862 x 10 -11 h' 3 + 1.25718 x 1(T 15 h 14 (38) 

The runway x-axis coordinate for the fast aircraft is then: 

Xfast ~ — COt y (hy a5t — hrcH ) (39) 

As with the TAS model, different values of h TCH , h FAF , and y can be used for the fast aircraft in Eq. (6), the 
Deas equivalent of Eq. (8), and Eq. (33) to Eq. (39) to accommodate runways with different threshold 
crossing heights, different elevations, or different glideslope angles. 

Finally, the required longitudinal separate at the FAF is: 

$FAF = %FAF ~ Xfast (40) 

If X FAF differs between the two runways, then use X FAF for the runway of the slow aircraft. Equation (40), 
like Eq. (15), assumes that the runway thresholds are not offset along the x-runway axis direction. To 
compute the separation compression use Eq. (16). 

Using the prior example of the slow aircraft with a final approach speed of 120 KT and a fast aircraft 
with a final approach speed of 130 KT, the results using the EAS model are: 
t sk)w = 148.1 seconds 

ti-fast = 135.9 seconds, (t sk)w - ti_ fast ) = 12.2 seconds > t dday so compute X D . fast 

Wei = 49.3 seconds 

D fast = 36699 ft 

h fast = 1948 ft 

X fast = -36222 ft 

^compress 2080 ft 

S faf = 2830 ft 


C. Comparison of Kinematic Models 

In comparison to the earlier example using the TAS model, the compression length increased from 2017 
ft to 2080 ft, a difference of 63 ft. So, how can the compression length be larger when the time of flight for 
both aircraft is shorter? The answer is that the fast aircraft maintains a higher altitude and, therefore, has a 
slightly higher true airspeed, even when the equivalent airspeed is equal to the slow aircraft. Depending on 
the initial separation, this effect can cause added compression that is greater than the reduced compression 
from the shorter flight time. It is one limitation of the EAS model. In the procedure, the fast aircraft uses a 
speed management algorithm that maintains separation through deceleration, effectively causing the fast 
aircraft to track the true airspeed of the slow aircraft. Therefore, until the fast aircraft decelerates to its final 
approach speed as an equivalent airspeed, the fast aircraft's equivalent airspeed should be slightly lower 
than the equivalent airspeed of the slow. One could partly compensate for this effect by giving the fast 
aircraft a slightly lower EAS for the constant segment. 

The above example places the fast aircraft about 148 ft above the slow aircraft when the slow aircraft 
reaches the FAF. If the slow aircraft has an EAS of 180 KT and the fast aircraft is matching the true 
airspeed of the slow aircraft, then the EAS of the fast aircraft would be 179.6 KT when the slow aircraft is 
at the FAF. Using this slightly lower EAS for V c for the fast aircraft would reduce the compression 
distance to 2042 ft, a change of 38 ft that leaves a difference between the two models of 25 ft. This 
adjustment does, however, require at least two iterations of Eq. (29) through Eq. (40). The first iteration 
approximates h fast in order to adjust V c for the fast aircraft using Eq. (21), and the second iteration 
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substitutes the adjusted V c . Otherwise, there is no simple mechanism to modify the EAS-based model to 
emulate true airspeed matching through deceleration. Nevertheless, the example results do demonstrate 
that the EAS-based model often provides the more conservative estimate of separation compression. 
Moreover, the EAS-based model remains sufficiently simple to program or to use in generating look-up 
tables; therefore, it can be used in real-time by flight deck or controller systems to generate custom values 
of the front-gate separation for specific aircraft pairs. 

The differences between the EAS kinematic model and the TAS kinematic model are very modest in 
the sea-level example. The difference in travel time does change substantially with runway elevation. 
Changing runway elevation to 1000 ft in the example above decreases t s i ow by 4.1 seconds compared to the 
TAS model, almost double the difference at sea level. However, the difference in separation compression 
is less sensitive to elevation since compression is a function of relative velocity between the aircraft pairs 
and the relative velocity changes more slowly with elevation given that each aircraft experiences nearly the 
same EAS-to-TAS conversion along the path. The difference in separation compression grows to only 70 
ft or, if setting V c of the fast aircraft to 179.6 KT to match the true airspeed of the slow aircraft at the FAF, 
the compression difference grows to 32 ft. 

The compression difference, however, is more sensitive to the difference in approach speed between the 
slow and fast aircraft. If the fast aircraft speed were increased to 140 KT on approach to the generic 
runway, then the difference in separation compression increases to 96 feet or, if setting V c of the fast 
aircraft to 179.6 KT, the difference increases to 65 ft. Thus, the EAS model becomes more relevant when 
analyzing the limit on the speed differential for a paired approach procedure. Also, as paired approaches 
are applied to runway separations as little as 700 feet, wake considerations drive the separation window to 
smaller lengths and require ever tighter performance margins. Ref. 4 estimated that the fast aircraft must be 
able to maintain position within a separation window of 490 feet to enable parallel, paired approaches 
without passing at 750 ft runway separation. The differences in separation compression estimates, 
therefore, become a greater percentage of the separation window as runway separation decreases and can 
influence the determination of required aircraft navigation performance and feasible aircraft pairings. 

D. Comparison of Kinematic Models to High Fidelity Simulation 

To validate the EAS model, its predictions of travel time from FAF to runway threshold were compared 
against a high fidelity simulation of a large civil transport. The high-fidelity simulation is one of the 
models available in the Cockpit Motion Facility at NASA Langley Research Center. 9 The high-fidelity 
simulation executed the approach using an autothrottle programmed to the speed schedule and an autopilot 
configured to follow localizer and glideslope signals. Furthermore, sensor errors were disabled so that the 
modeled aircraft would better conform to the approach path. However, conforming to the speed schedule 
was more challenging. The high fidelity simulation doesn’t have an autothrottle mode to mimic the 
constant deceleration in the kinematic models. Moreover, deceleration depends in part on landing gear 
deployment, flaps settings, and speed brake settings, all of which require pilot action. Therefore, the 
simulation uses a pseudo-pilot algorithm to operate the levers for gear, flaps, and speed brake. In this 
simulation, the speed-brake deployment is a one-time action at deceleration start, and the lever position is a 
function of the the target deceleration. The speed brake is then retracted prior to achieving the final 
approach speed. Table 2 shows the predictions of travel time using the TAS model, the EAS model, and the 
simulation. The approach was to KSFO RWY28L (TCH = 57 ft, runway elevation = 13 ft): 


Table 2. Comparison of Approach Times. 


Approach 
Speed (KT) 

Time from FAF to TCH (seconds) 

TAS 

model 

EAS Model 

High 

Fidelity 

Simulation 

120 

149.34 

147.58 

146.66 

130 

140.55 

138.87 

138.58 

140 

132.86 

131.25 

131.34 

150 

126.06 

124.52 

124.84 

160 

120.00 

118.52 

118.12 
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Except for the 120 KT approach speed, the estimates of travel time from the FAF to the runway 
threshold are nearly the same (within + 0.4 seconds) for the EAS model and the high fidelity simulation. 
For the 120 KT approach speed, the high fidelity simulation arrives 0.92 seconds earlier than estimated by 
the EAS model. The 120 KT approach speed requires a deceleration near the limit of what the modeled 
aircraft can achieve when autothrottle is engaged, and the actual speed of the aircraft in the high fidelity 
simulation is closer to 121 KT at the SAP. When the EAS kinematic model was re-evaluated with an 
approach speed of 121 KT, then the estimated time to threshold changed to 146.47 seconds, which is a 
closer match to the high fidelity simulation. 


E. Variations to the EAS Kinematic Model 

1. Deviations from Standard Day Density 

The EAS model, as presented in section III.B, is tailored to the standard day profile for density in the 
1976 Atmosphere Model. Actual density that departs from the standard day profile will also adjust the 
compression distance. One can employ the methods used to generate Eq. (24) against adjustments to the 
standard day. For example, when adjusting density by changing the sea-level temperature, the coefficients 
of Eq. (24) can be formulated as functions of the sea-level temperature (T 0 ) in degrees Rankine: 


DeasQ 1 , T 0 ) 


0.00379429 _ 
cscy Qi h z 


(5.08743 x 10" 6 + 1.21288 x 10" 11 T 0 ) „ 
+ - 7*2 -h 3 ) 

‘0 


(41) 


Unfortunately, the conversion back to geometric altitude, i.e. the real root of Eq. (41), is a long formula: 5 
h = (0.00126476 T 0 + ((-3 .63534x1 O' 7 - 1.52814xl0 10 T 0 ) T 0 2 ) / (6.98812xl0 10 h' T 0 2 - 
6.44784xl0" 8 T 0 3 + 3.33204xl0 14 h' T 0 3 - 4.14183xl0 12 T 0 4 + 3.97192xl0 19 h' T 0 4 + (4 
(8.6561xl0' 7 + 3.63865xl0 10 T 0 ) 3 T 0 6 + T 0 4 ((-6.44784xl0 -8 - 4.14183xl0 12 T 0 ) T 0 + h' 

(6.98812xlO' 10 + (3.33204xl0‘ 14 + 3.97192xl0' 19 T 0 ) T 0 )) 2 ) m ) 1/3 + 0.264567 (6.98812xlO' 10 h' 

(42) 

T 0 2 - 6.44784xl0' 8 T 0 2 + 3.33204xl0' 14 h' T 0 3 - 4.14183xl0' 12 T 0 4 + 3.97192xl0' 19 h' T 0 4 + (4 
(8.6561xl0‘ 7 + 3.63865xl0' 10 T 0 f T 0 6 + T 0 4 ((-6.44784xl0‘ 8 - 4.14183xl0' 12 r 0 ) T 0 + h' 

(6.98812xl0 10 + (3.33204xl0 14 + 3.97192xl0‘ 19 T 0 ) T 0 )) 2 ) 1/2 ) 1/3 ) / (5.08743xl0' 6 + 

1.21288xlO" 10 T 0 ) 

While leaving T 0 as an input parameter, Equation (42) cannot be simplified further using a truncated Taylor 
series in h’ [as was done to obtain Eq. (38)]. It may be possible to eliminate some of the terms as negligible 
but this exercise is left to the reader. Despite its length, Eq. (42) remains a workable formula for coding 
and executing on computer. Alternatively, one could apply a series of sea-level temperatures to Eq. (41) 
and Eq. (42) in order to produce a series of functions in h' alone. Furthermore, the resulting series of 
inverse functions [i.e. those derived from Eq.(41)] could be further simplified using a truncated Taylor 
series. The resulting coefficients for the series of formulas could then be stored in a look-up table. This 
alternative would require significantly less computation at the cost of additional memory. 

2. Adding Winds 

The model presented also does not include winds. Winds speeds are true speeds that are often modeled 
as functions of altitude, i.e. W(h). Thus, in the presences of winds, Eq. (23) changes to: 

/ KT \ dr ( KT \ dh 

= (j^h +wm = = v ‘-f m <43) 


§ As solved and simplified by Mathematica™ Release 8.0 
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S i,=,= { cscr ?h)S 


fpsJJ (V EAS f(h)-W(h)) 


dh 


(44) 


The integral for dh can be estimated using a Taylor series expansion as was done for Eq. (24). However, a 
similar algebraic formulation cannot be derived for the deceleration segment because the time and altitude 
terms cannot be segregated into separate integrals: 


V, 


TAS 


{ 


KT \ dh 


CSCY fps)'dt +W ^ = ^ EAS + Ueas ^ ^ h>> 


(45) 


(We as + a EAS t ) f(h ) - W(h)) dt = (esc Yj^j 


dh 


(46) 


If one makes a simplifying assumption that the wind speed is an equivalent speed rather than a true speed, 
Eq. (44) and Eq. (46) change to: 


( KT \ f 1 

f V CSCK /WJ (v EAS — w (h))/ (/i ) dh 


(47) 


n ( KT\ 

( Veas + a EAs t - w (h))f(h) dt = (esc y — j dh (48) 

Again, the time and altitude terms for the deceleration segment cannot be segregated into separate integrals 
except in the case where the wind speed is constant with altitude. In this case, the resulting equations are 
identical to Eq. (24) to Eq. (34) with (Veas - W) substituted for Veas- Therefore, this is the only option for 
generating an algebraic formulation that includes winds. 

3. Adapting the Model to Offset Runways 

As previously discussed, formulation of the models as functions of altitude and glideslope angle 
facilitate the model’s adaptation to runway pairs whose approach geometries may differ by threshold 
crossing height, runway elevation, final approach fix altitude, and glideslope angle. Accommodating 
runways that are offset in the longitudinal direction requires further adjustment because the equations for 
the procedure termination altitudes [Eq. (4) and Eq. (6)] implicitly assume no offset in the x-axis origin 
between runways. These equations derive from the assumption that the procedure ends when an aircraft 
crosses the threshold and that longitudinal location is the same whether the slow or fast aircraft crosses the 
threshold first. With offset runways, the procedure termination altitudes, h PTS and h PTF , now depend on 
which aircraft is expected to cross the threshold first. The paired-approach without passing is the simpler 
case because the slow aircraft is always expected to cross the threshold first. For this case, h PTS remains 
unchanged (using Eq. (4) where D ws is zero). However, Eq. (6) must be modified to include the runway 
offset distance (D offset ) in determining h PTF : 

h PTF = h TCH + (D cs + D off se t ) tan y (49) 

However, in the case that allows passing, the fast aircraft can cross the threshold first only if the runway 
offset distance (D offset ) is shorter than wake-safe distance (D ws ). Otherwise, the slow aircraft must cross the 
threshold first if it is to avoid a wake encounter after the first aircraft passes. In this case, the fast aircraft, 
though ahead of the slow aircraft, is also behind its runway threshold when the procedure terminates. In 
the case where the fast aircraft reaches the threshold first, h PTF remains unchanged (using Eq. (6) with D cs 
equal to zero) and Eq. (4) must be modified to account for D offset : 

hpTS = h T CH + (Avs — ^offset) tan Y (50) 

In the case where the slow aircraft must cross the threshold first, then h PTS remains unchanged (using Eq. 
(4) where D ws is zero) and Eq. (6) is adjusted for the runway offset as follows: 
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h PTF — h TCH + ( D 0 ff Set — D ws ) tan y 


(51) 


The equations that compute the runway x-axis distance of the fast aircraft, required separation at the 
FAF (S FAF ), and longitudinal compression (d compress ) also need adjustments for offset runways. Since the 
separation of the fast aircraft relative to the slow is the quantity of interest, the recommended approach is to 
adjust the equations to use the runway coordinates for the slow aircraft as the reference. For all scenarios, 
D offset is added to Eq. (10) in the TAS model and Eq. (40) in the EAS model. In the case of the paired 
approach without passing, no other adjustments are necessary. However, the TAS model requires 
additional adjustments for the paired approach with passing that depend on which aircraft is expected to 
cross the threshold first. If the fast aircraft will cross the threshold first, then the TAS model must also add 
D offset to Eq. (14). If the slow aircraft will land first, then the TAS model adds the wake-safe distance D ws 
instead of D 0 ff se t to Eq. (14). 

F. Model Limitations in Addressing Aircraft Performance Variability 

The equations in this paper model ideal performance; however, the actual performance of aircraft 
varies. First, the model does not account for path deviation (lateral or along-path flight technical error) 
which can modify the flown distance (and therefore flight time) from FAF to touchdown. The model also 
does not directly address the variation in the actual altitude where the slow aircraft initiates deceleration. 
Late deceleration will shorten the compression distance; early deceleration will lengthen it. Nevertheless, 
one can use the model to examine the bounds of late or early deceleration by assigning the altitude of the 
late or early deceleration to the variable h FAF for the slow aircraft. Likewise, actual deceleration not only 
can vary from the ideal, but it isn’t a square impulse; it ramps up and ramps down. As demonstrated by the 
comparison with a high-fidelity simulation in section III.D, variability of deceleration appears to have small 
effect. Even so, one can reformulate the equations to accept a user-defined constant acceleration as an 
input in order to explore the variation bounds for deceleration. This would require additional equations to 
determine how much of the true (TAS model) or D EAS (EAS model) distance from the FAF to the 
procedure termination height (PTS or PTF) is covered by the deceleration segment and how much remains 
to be flown under the final approach speed. These equations are necessary to compute both t s i ow and t Ffast . 
The equations for X D _ fast (TAS model) and D fast (EAS model) are already dependent on an assigned 
acceleration (that of the slow aircraft) and would not need modification. For example, to compute t s i ow with 
a user defined acceleration, Eq. (2), Eq. (4), and Eq. (5) would be replaced with 


_ AK _{V f -V c ] 

) 

(52) 

tdecel 

U slow U slow 


(v f + v c ) 

n — t v — t v J 

u decel L decel v avg L decel ^ 

_ (v? ~ v c z ) 

^Uslow 

(53) 

(DpAF^PTS ~ ldecel ) ( D(h FAF ,y ) — 

D(h PTS ,y ) — D dece i ) 

(54) 

v f 

Vf 

tsiow ~ tdecel f/inai 


(55) 


where D(h,y) is the function for the TAS model [see Eq. (1)] or the EAS model [see Eq. (24)] that converts 
altitude and glideslope angle to glidepath distance from seal level. 

Because much of the distance from FAF to threshold is traveled independently, deviations in actual 
aircraft speed from the scheduled speed can have a large impact on the decrease in longitudinal separation 
between the aircraft. 3 The TAS and EAS models are both derived using constant speeds and the EAS 
model, in particular, cannot be modified to inject time or distance dependent variations in speed during a 
run. However, one can assume that, over the course of the approach, the effects of the random component 
of speed variation largely cancel out and that only the bias affects outcome. Thus, one can explore the 
bounds on speed deviation by applying constant biases to V f and V c for the slow and fast aircraft. The 
worst-case decrease in separation occurs when the negative limit on the velocity bias is applied to the slow 
aircraft and the positive limit on velocity bias is applied to the fast aircraft. 4 Reference 3 provides an 
example of applying velocity biases to the TAS model to determine the required separation at the FAF that 
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is robust to velocity bias. Similarly, Ref. 4 provides an example of applying velocity biases to the EAS 
model to determine required separation at the FAF that is robust to velocity bias. 

IV. Conclusions 

A set of algebraic formulas were derived to estimate the travel time and separation loss of aircraft 
participating in a paired approach when the speed schedule of those aircraft is specified using equivalent 
airspeeds. This model was validated against a high-fidelity model of a civil transport flying the same path 
under autopilot and autothrottle control. Estimates of time to travel from FAF to threshold generally 
matched the high fidelity simulation within + 0.4 seconds. This paper also presented extensions of the EAS 
kinematic model to incorporate differences in sea-level atmospheric density from standard day (using the 
sea-level atmospheric temperature as a model input) and to incorporate constant winds though this 
extension requires treating the wind speed as an equivalent speed rather than a true speed. However, the 
EAS model could not be extended, algebraically, to wind speed as a function of altitude. The model also 
assumes an ideal path and does not address the path deviations caused by flight technical error. 
Nevertheless, the model inputs or its equations can be adjusted to explore the deviation bounds in the 
altitude of deceleration start by the slow aircraft, deviation bounds on the average deceleration of each 
aircraft, and bounds on bias in following the speed schedule. 

When the EAS model is compared to similar algebraic formulas derived from the assumption that speed 
schedule is flown as true airspeed (TAS model), the differences in estimates of travel time and separation 
loss are modest. In the examples presented here, travel time differs up to 4 seconds, and separation loss 
differs by up to a 100 feet. Nevertheless, the EAS model results are more conservative than the TAS 
model. Moreover, paired approaches demand ever tighter performance margins as they are applied to 
runways separated by as little as 700 feet. Under such performance constraints, these modest differences 
have the potential to be significant determinants of required aircraft navigation performance and feasible 
aircraft pairings. The equations also remain simple enough to implement in real-time decision support 
tools as a higher fidelity alternative to the TAS model. 

References 

1 Johnson, S. C., Abbott, T. S. , Guerreiro, N. M. G, Lohr, G. W., Volk, P., and McKissick, B. T., 
“Simplified Aircraft-Based Paired Approach: Concept Definition and Initial Analysis,” National 
Aeronautics and Space Administration, NASA/TP-2013-217994, 2013. 

2 Bone, R., O. Olmos, and A. Mundra, “Paired Approach Operational Concept,” Proceedings of the 20th 
Digital Avionics Systems Conference, Daytona Beach, FL, 20001. 

3 Perry, R. B., Madden, M. M., Torres-Pomales, W. and Butler, R. W., “Analysis and Simulation of the 
Simplified Aircraft-Based Paired Approach Concept With the ALAS Alerting Algorithm in Conjunction 
With Echelon and Offset Strategies,” NASA/TM-2014-218151, 2014, publication pending. 

4 Torres-Pomales, W., Madden, M. M., Butler, R. W. and Perry, R. B., “The Simplified Aircraft-Based 
Paired Approach with the ALAS Alerting Algorithm,” NASA/TM-20 13-2 17804, 2013. 

5 Hammer, J., “Case study of paired approach procedure to closely spaced parallel runways,” Air 
Traffic Control Quarterly, Vol. 8(3), 223-252, 1999. 

6 Eftekari, R. R., Hammer, J. B., Havens, D. A., and Mundra, A.D., "Feasibility analyses for paired 
approach procedures for closely spaced parallel runways," Integrated Communications, Navigation and 
Surveillance Conference (ICNS), 2011 , pp. 15-1 - 15-14, 10-12 May 2011. 

7 Federal Aviation Administration, "Automatic Dependent Surveillance — Broadcast (ADS-B) Out 
Performance Requirements To Support Air Traffic Control (ATC) Service; Final Rule," Federal Register 
(National Archive and Records Administration) 75, no. 103, May 2010, pp. 30160 - 30195. 

8 National Oceanic and Atmospheric Administration, National Aeronautics And Space Administration, 
and United States Air Force, "U.S. Standard Atmosphere, 1976," NOAA-S/T 76-1562, U.S. Government 
Printing Office, Washington, D. C., 1976. 

9 Chung, V. I., and Hutchinson, B. K., “A Unique Software System For Simulation-To-Flight 
Research,” AIAA Modeling and Simulation Technologies Conference & Exhibit, Montreal, Canada, AIAA 
2001-4057, August 6-9, 2001 


16 



